Thawing f{R) cosmology 
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We consider Brans-Dicke (BD) scalar tensor theory in the conformally transformed Einstein frame. 
In this frame BD theory behaves hke an interacting quintessence model. We find the necessary 
conditions on the form of the potential V{ifi) in order to have thawing behavior. Finally, by setting 
the BD coupling constant u = 0, the metric f{R) gravity has been considered in the Einstein frame. 
Assuming the existence of thawing solution, some necessary conditions for f{R) gravity models have 
been derived. 
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INTRODUCTION 

One of the proposals for explaining the present accel- 
erated expansion of the universe [H is modifying Ein- 
stein's theory of gravity by introducing corrections to 
the Einstein-Hilbert lagrangian. These theories, called 
"modified gravity theories" follow this idea that the 
accelerated expansion of the universe may be has a geo- 
metric interpretation instead of adding the exotic forms 
of energy sources, dubbed "dark energy" (sj. In the other 
words, in this perspective the dark energy is a manifesta- 
tion of a modified gravitational interaction rather than a 
new form of energy density. The situation is reminiscent 
of the problem of precession of Mercury's orbit. In the 
mid-nineteenth century the anomalous behavior of Mer- 
cury firstly was attributed to some unobserved (" dark" ) 
planet in the solar system while it was mainly due to 
the failure of Newton's theory of gravity in the strong 
gravitational field regime. In this view, it seems that 
as long as the dark energy particles [4,] have not been 
observed directly, the "geometric" candidates have im- 
portant role. The simplest form of the modified gravity 
theories can be obtained by replacing the Ricci scalar R 
with an arbitrary general function f{R) in the Einstein- 
Hilbert action, usually called f{R) theory of gravity. For 
a recent review of this theory see . 

Metric f{R) gravity model is dynamically equivalent 
to a BD scalar tensor theory with coupling constant 
a; = By using this equivalence, one can eas- 

ily find the prediction of metric f{R) gravity for the 
PPN parameter "fppN- This parameter in the BD the- 
ory has the form ^ppn = (1 -I- w)/(2 -I- w). Thus the 
value of this parameter in the metric f{R) gravity is 
1/2 which it is not in agreement with the experimen- 
tal bound \^ppN - 1| < 2.3 : 10"^[7j. However, con- 
sidering this model in the Einstein frame has some sat- 
isfactory features. For example, f{R) gravity can dis- 
play the chameleon behavior in this frame which helps 
to relax the weak field limit problem of f{R) gravity [8|. 
Chameleon effect is firstly interpreted using the scalar 
tensor framework of dark energy In this theory the 
effective mass of the scalar field is a function of the curva- 
ture of space-time and consequently it can be large at the 
solar system and small on the cosmological scales. This 



behavior appears in the minimally coupled scalar ten- 
sor theory if there exists an energy transfer between the 
dark energy fluid and the ordinary matter fluid. Since the 
quintessence model |lol | is a minimally coupled scalar ten- 
sor theory, the chameleon mechanism can be appeared. 
On the other hand, metric f{R) gravity theories are con- 
formally equivalent to models of quintessence in which 
matter is coupled to the dark energy, thus the chameleon 
effect can occur in the conformal frame [11 . 

The noninteracting quintessence models can be di- 



vided into two categories [Uj. "Freezing" models: in 



these models the equation of state parameter of dark en- 
ergy, ujip^ has an arbitrary value initially and decreases 
with time and asymptotically approaches —1. "Thaw- 
ing" models: these models have a value oi lu^ ^ — 1 
initially, and it increases with time. There is a subset 
of freezing models which display tracking behavior [l^ . 
In the tracking models, lo^ has an arbitrary value ini- 
tially and it is nearly constant during the tracking era. 
When the tracking era terminates then oj^ decreases and 
asymptotically approaches -1. The important feature of 
these models is that the evolution of the scalar field is 
insensitive to the initial conditions and the dark energy 
density drops with a slower rate than the matter energy 
density and finally overtakes it. Albeit, these models can 
not provide a solution to the so-called coincidence prob- 
lem because other fine-tunings are needed on the free 
parameters of these models in order to have an appropri- 
ate amount of dark energy compatible with observation 
in the present davs[l3j. 



In our recent paper 1J| we derived some conditions 
for existing the stable tracker solutions in the Einstein 
frame of metric f{R) gravity models. It is found that 
the tracker solutions with —0.361 < uj^ < 1 exist if 
< r < 0.217 and f^\nf\R) > 0, where F is a di- 
mensionless function defined by relation pT|) in the next 
section. The main purpose of this paper is to find out the 
necessary conditions for the existence of thawing behav- 
ior in the Einstein frame of metric f{R) gravity theories. 

The outline of this paper is as follows: In section II 
we start with BD scalar tensor theory (with an arbitrary 
uj). As mentioned before, this theory behaves as an inter- 
acting quintessence model in this frame. We derive some 
necessary conditions on the form of the potential V{ip) in 
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order to lead to the thawing behavior for uj^. In section 
III, by setting a; = in the results, we present a gen- 
eral description of the behavior of the thawing f{R) in 
the Einstein frame and finally conditions on the form of 
f{R) gravity have been derived. Throughout this work 
we have chosen the unit SttG = c = 1, the metric sig- 
nature is (H ) and the universe is assumed to be 

spatially flat. 



THAWING NONMINIMAL QUINTESSENCE 

The effective action for BD scalar tensor theory is given 

by 

(1) 

where R is the Ricci scalar, [/(<&) is the potential of the 
scalar field and Cm{gnu) represents the matter lagrangian 
density. Note that all tilded quantities are in the Jordan 
frame. The coupling constant lo should be large to pass 
the experimental testes. The observational constraint on 
w is |a;| > 40000 [tI. Under the conformal transformation 

QtLv = e^'-^^v (2) 

where In $ = Qip and C — ^rrhr- , one can obtain the 
Einstein frame action 



3+2w ' 



Sr — 



J ^ d^x [R - ^r^^^ 



(3) 



where V[ip) = e~'^^'^U{<^{(p)). We see that in the Ein- 
stein frame the scalar field couples conformally to mat- 
ter via the function e"'''^ but couples minimally to the 
gravity sector. For a spatially flat FRW universe, the 
modified Friedmann equations are given by 



1, 



+ Pm) 



H = - 2 + ^m)Pm + (1 + '^ip)Pip\ 

and the equation of motion of the scalar field is 



/3(1 - iuJm)p 



(4) 



(5) 



where = y f C> is the equation of state parameter 
of the ordinary matter with the energy density p„i in 
the Einstein frame. Also = \ip^ + V{ip) and — 



i,;,2 



^ip~ — V{}p) represent the energy density and pressure of 
the dark energy respectively. The conservation equations 
of the scalar field fiuid and the cosmic fiuid are 



Pm + iH{l + UJm)Pra = - V g/3<y5(l - SWm)^^ 



(6) 



and the energy density of matter pm , pressure Pm , cosmic 
time t and the scale factor a are related to their Jordan 
frame counterparts through [l^ 

Pm = e.^'^'^'^PmiPm = e~^'^'^p„ , = e^di,a = e^d (7) 

During the matter dominated era, by using equation ([5]), 
one can introduce an effective potential as follows 



Veff{ip) = Viip)+p*e"Vl0'P 



(8) 



Where p* is a conserved quantity in the Einstein frame 
0, which is related to pm via the relation pm = 

Since the late time evolution of the universe is of inter- 
est here and also our main purpose is to explore the role 
of the interaction term (which is nonzero for the mat- 
ter component), we neglect the radiation component and 
assume that the universe contains only dust and dark en- 
ergy. It is interesting to note that the interaction term 
is commonly assumed to be zero in the radiation dom- 
inated era, but recently Cembranos and et al (l6| have 
shown that this interaction term can lead to strong im- 
pact on cosmology in the radiation dominated era due 
to the finite temperature radiative corrections. In the 
other words, there exists another source term for scalar 
field given by the conformal anomaly which leads to a 
nonzero trace of energy momentum tensor in the radi- 
ation dominated era (note that the RHS of ([5]) is the 
trace of energy momentum tensor) . Considering the con- 
formally coupled scalar field with a quadratic coupling 
function and vanishing potential, the above effect leads 
to a temporary contracting phase in which the temper- 
ature increases [lij. However, as mentioned before, we 
aim to study here the late time evolution of the universe 
and so we assume that the universe is filled with non- 
relativistic matter. 

Following reference [tH ]. we introduce the variables x, 
y and A defined by 



V6' 



y 



X 



V 



(9) 



where the prime denotes the derivative with respect to 
In a. By these definitions, it is an easy job to show that 
the equations ^ and ([5]) become 



x' = -3x + \y' + \x{\ ^x"- y") + /?(1 - x" - y") 



y' = ^-^y^^y + + ^'^ ^ 2/^) 
X' ^ -Vex^ir -i)x 

where 

dip dip 



(10) 



(11) 
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For thawing models 



-1 and so 7 = 1 + <C 1. 



Thus it is convenient to express the above equations with 
respect to 7 in order to exploit its smallness by expanding 
quantities to the lowest order in 7. Also we assume that 
(/? > (x' > 0). This assumption, as considered in [14| . 
is necessary to have an increasing dark energy density 



For uncoupled quintessence, where f3 is zero, the RHS 
of equation p7)) is approximately constant and moreover 
it has an small amount (note that 7 <C 1) for thawing 
solutions . Thus if (■^)^ ^ 1 then the thawing behav- 
ior can occur [l^. In the general case where /3 is not 
zero, then the RHS can not be regarded as a constant. 



parameter i.e. fl^, > 0. However, the results can be Since fl^ is appeared in the denominator, hence the sec- 



generalized to the opposite case (x' < 0). Now, by using 
y'^ and 7 = 2a;^/i7^, one can rewrite the 



equations ()10|1 in terms of 7 and f2(^ as 



7' = (2-7) -37 + ^^37^^ 




r!^) (12) 



n'^ = 3{i-np)(^{i~^)n^ + ^y/2^^ (13) 



A' 



-V3A2(r- 1)^^ 



(14) 



It is clear from equation (fT^]) that for thawing models 
il'^p during the cosmological history of the universe 
(0 < flip < 1). Thus by using (|13|1 we can write equation 
(fT2l) as follows 



d7 _ (2-7)(-37 + Av/3^+y^/3(l-17^)) 



3i)^(i-i7^)(i-7 + fy^) 



(15) 



This equation is obtained earlier in [IS"! in which, the 
non-minimal quintessence with nearly flat potentials has 
been considered. Equation (|15p is not a simple differen- 
tial equation and for solving it, we will make some as- 
sumptions which are satisfied for thawing models. First 
assume that 7 <C 1, by retaining terms up to the first 
order in 7, the equation (|15p takes the following form 



c?7 
dD,,„ 



-27 




A7 



2AV7 



y3a;(i - n^) 



2V 



(16) 



3/2 



Another useful equation can be obtained by using the 
equation of motion of the scalar field ^ 



ond term in the RHS is dominated initially and has a 
large value. Thus the LHS can not have a small value 
as well as can not be a constant. When gets larger, 
the effect of the interaction becomes weaker. Thus, at 
late times, the nearly flat region of the potential leads 
to the thawing uncoupled quintessence. So, unlike the 
noninteracting quintessence model, nearly flat potentials 
can not lead to the thawing behavior when an explicit en- 
ergy transfer between the scalar field fluid and the matter 
fluid exists. In this mentioned in |18l | , with nearly 

flat potentials, firstly increases with time and then, 
when the interaction becomes weaker, it decreases and 
approaches asymptotically to a value near —1. The be- 
havior of Lo^p with nearly flat potentials has been plotted 
in Fig[T]by solving equation (fT5|) numerically. Note that 
this behavior is due to the special form of interaction 
which appeared here (i.e. tppm)- 

Now we are ready to make the second assumption. 
Taking into account equation (jlSp and assuming that the 
value of the term within the bracket to be approximately 
constant for thawing solutions, this equation gives 



A~ An- 



3 



Ao 



(19) 



where Ao is a positive constant. Hereafter we shall refer to 
this equation as the "thawing condition". For potentials 
in which — ^ decreases as ip increases (F > 1), the LHS 
of the equation is increasing. On the other hand, 
the RHS is increasing because il^p increases. Hence, the 
thawing condition can not be satisfied. Thus, the thaw- 
ing condition shows that it is necessary A increases with 
time when tp and are increasing, i.e. 



F < 1 



(20) 



It is clear from the thawing condition that if /3 = then A 
is nearly constant and so F ~ 1. One can find other sim- 
ple conditions on the form ofV{ip) by using the thawing 
condition. For this purpose, let us rewrite (fT9|) as follows 




7' 



^ 1 - r!., 



37(2-7)^ 2-7 



(17) 



This equation can be written to the first order in 7 as 
follows 



(18) 



3 



(21) 



It is clear from this equation that the second term in the 
RHS is dominated initially and so ^ > 0. As men- 
tioned before, the interaction becomes weaker at late 
times. Thus, there exists a time i*, at which 



Ao + 



/3 



/3 



(22) 



¥>/ t=t* 
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FIG. 1: Numerical solutions of ([T5|) 
when (3 — 0.5 for (top to bottom) A 



for nearly flat potentials 
= 1, A = 0.5, A = 0.1 and 



A = 0.01. Assume that lj^ 



-1 at = 0.001. 



FIG. 2: Our analytical result for lj^ as a function of for 
[3 = 0.5 and Q^o = 0.7 for (top to bottom)Ao = 1, Ao = 0.9 
and Ao = 0.8. Also, as an initial value, it has been assumed 
that at Q.^ — 0.001 7 is zero. (Sl^o is the current value of n,p) 



At this time 



< 0. Thus, it is 



f ^ and after it ^ ^ 
necessary that the potential has a maximum in order to 
have thawing behavior. In the other words, a value of 
the scalar field, lo* should exist such that 



=0 



dv 



(23) 



dip 



2 \ v* 



< 



Now let us to justify the thawing condition. By these 
assumptions (7 <C 1 and (fT9|l ). equation (fTH]) takes the 
very simple following form 



c?7 



2Aj 



2A0V7 







(24) 



in which A= 1 + f3y ^Xq. Note that if (3 — 0, then this 
equation becomes precisely the equation obtained earlier 
by Scherrer and Sen [l9j . This differential equation has 
an exact solution as follows 



2A 



\/3(l + 2 A) 



(25) 



AA 



A,- + A,n^)f 



where xo is an integration constant depending on the 
initial conditions and 2F1 is the Gauss Hypergeometric 
function. This equation gives an analytical expression 
for the state parameter of dark energy as a function 
of its density parameter for the thawing non-minimal 
quintessence model. It generalizes the result obtained 
in [19] for thawing minimally coupled scalar field. The 
behavior of uj^p as a function of Sl^ has been shown in 
Figl2] for various values of Aq. The /3 has been chosen 
to be 0.5, the value will be used in the next section in 
the case of f{R) gravity models. In fact, the value of 
Aq should be such that lo^ has a value near —1 today. 
For confronting the model with observational data, it is 




FIG. 3: The solid curve represent the exact solution of (|13p 
when 7 = 0, i.e. the equation (|26|l . assuming fi^o = 0.7. 
The dot dashed curve is the numerical solution of (|13|l with 
Ao = 2, the dotted curve is for Ao = 1 and the dashed curve 
is for Ao = 0.8. 



needed to express 7 and fi^ in terms of cosmic red shift 
or cosmic scale factor. By substituting the solution P5|) 
in the equation (jlSp , we obtain a differential equation 
for U,^ in which the Gauss Hypergeometric function is 
appeared. Here, we have solved it numerically and the 
result has been compared with the exact solution of the 
equation (|13p when 7 = 0, in Figl21 Thus as it is clear 
from Fig[3l the difference between these solutions is small 
when Aq '--^ 1 and consequently one can use the solution 
of (fT5|) when 7 = i.e. 



(26) 



in order to find out an approximated expression for 7 as 
a function of a. 



THAWING f{R) 



Now let us consider f{R) gravity in the Einstein frame. 
It is sufficient to set w = (and so C, — -i /|) in equation 



5 



(|25)) in order to have thawing behavior in the Einstein 
frame. Also in the context of f{R) gravity, the scalar 
field ip is related to the curvature scalar of the Jordan 
frame as follows 



if : 



/3721n/^(i?) 



V{ip) = (i?/^ - /)/2/| 



(27) 



R dR 

chosen near to 1 in order to have uj, 



where /s = As it is clear from Fig[21 Aq has been 

^ art 

near —0.9. In this 



case, as it has been shown in FiglJl fl^p evolves as the 
dark energy density parameter of ACDM model in which 
uj^ is always approximated to —1. 

Now, we want to find out explicit conditions on the 
form of the function /(i?) in order to have thawing be- 
havior in the Einstein frame. Taking into account equa- 
tions (PO)) and (P7)) . one can easily verify that if /^^^ > 
then 



Rfl - Wrr + fR)f < 



and for < 



Rfl - Wrr + fk)f > 



(28) 



(29) 



Also, by using and ^7} . it is necessary that the form 
of f{R) be such that there exists R* for which 



^ - -tIr' 

JR 

fk I _ 



(30) 



R* > 



RR 



Note that for having nonsingular conformal transforma- 
tion (equations ([77]) and ^) we have assumed /j| > 0. 
Equations ([28l) - ((30)) are the necessary conditions on the 
form of f{R) for raising to the thawing behavior and they 
are not sufficient conditions. 

Now, let us to find out an explicit example for thawing 
potentials. For this purpose, assume that 



dip 



where, a is a positive constant. This assumption leads 
to the following form of dark energy density parameter 



fl^p = 1 — e 



(32) 



which is an increasing function of (p and ijj is an integra- 
tion constant. Using cp' = ^yS^yfl^ and equations ([13]), 
(ISTIl and 111), we obtain 



(33) 



where B = 3- 2yf2a(i + 2a^f?. If a is a small quantity 
(a < 1), then the solution with minus sign can yield 



-0.85 



-0.90 



-0.95 



-1.00 



0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 

FIG. 4: ojip as a function of for fi — 0.5 for (top to 
bottom)Q = 2/7, Q = 1/4 and a = 2/9. 
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FIG. 5: Region of parameter space compatible with the obser- 
vational constraints — 1 < uj^q < —0.9 and 0.6 < Q^o < 0.8 
for P = 0.5. 



to the thawing behavior. For seeing this, uj^p has been 
plotted in FiglUfor various values of a. By substituting 
equation (|33p into equation ^T7\ and expanding the RHS 
of equation ^T7\ to the second order in a (note that we 
have assumed that a is small), we get 



1 dV 

Vdn^ 



1 +6»y 

2(1 -r!^) 



+ Oia^) (34) 



(31) which has the following solution 



|q/3 



(35) 



where Vq is a positive integration constant. It is obvious 
from this that V has a maximum and so it is consistent 
with our pervious results. By setting -0 to be zero and 
using equation (|32[) . let us rewrite equation (j35[) as a 
function of the scalar field as follows 



(1 



n 



(36) 



This potential satisfies the conditions ([20]) and ([23|) and 
it is a two parameter potential (Vq and a). The parame- 
ter a should be small and for —1 < uj^^q < —0.9 it should 
be < a < 0.23. Thus, the only free parameter in this 
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model is Vq. As mentioned before, this free parameter 
should be fin-tuned by using the observational data. The 
observational fact is that the energy density of dark en- 
ergy and the energy density of cosmic matter fluid are 
approximately in the same order. Since the potential has 
been obtained with respect to il^ (equation ([55]) ) it is 
easy to make an estimation on the values of Vq to repro- 
duce the acceleration expansion. Albeit, we assume that 
the major contribution to the energy of the scalar field is 
due to the potential term (note that this is the case for 
all thawing potentials). The density parameter of dark 
energy is 



V 
3IP 



By using this equation and ((35| we obtain 



Vo 
ml 



I a/3 



(37) 



(38) 



where Hq is the current value of the Hubble parameter. 
By taking into account that the current value of n^po 
satisfies the bound 0.6 < fl^o < 0.8, we have plotted 
the region of parameter space able to cover the above 
observational constraints, in Fig. [S] This region varies 
from a = to a = 0.23 and from Vq ~ 0.6 to Vq « 0.95, 
where Vq is a dimensionless variable defined as follows 



^^0 = 



Vq 
3Hl 



(39) 



Note that for making a more precise estimation one 
should use numerical solutions of the field equations and 
taking into account the effect of the kinetic term of en- 
ergy density of the dark energy, see the third paper of 
[ist and also for more details. 

Now, for finding the corresponding f{R) function, as- 
sume that /(i?) differs from Einstein's general relativity 
by a small perturbation as follows 



/(i?) = i? + e*(i?) 



(40) 



where e is a very small parameter. By substituting this in 
equation (P7)) . using ([55)1 and taking (3 = 0.5, one reaches 
to the following first order differential equation up to the 
first order in e 



which has the solution of the form 



where 



a — V6 
^ = 2K)(1-^) 



12 



(41) 



(42) 



(43) 



375 



Since a is a positive constant, n is a negative real number. 
Thus, the perturbation procedure is valid if the curvature 



< R. 



of space time is sufficiently high such that 

Consequently, the model can lead to the thawing 
behavior only in the beginning of matter dominated era. 
However, for larger values of tp, let us expand the poten- 
tial ([36]) again. Before proceeding, we expect that vE'(i?) 
contains some terms of R with powers smaller than n. 
Because such a term can have effect in the late times 
(large ip), where the curvature is small, while it can be ne- 
glected compared with i?" where the curvature is larger. 

If If is large enough such that e~° ^ 1 then one can 
write the potential ([36|) as follows 



V^Voe 



3a 



(44) 



It is easy to show that the f{R) function corresponding 
to this potential is 



f{R) ~ i^R'' 



(45) 



where 



7V6a~8 



a 



2Vbe3" (1 



(46) 



(2yoe3«'(2 



It is possible to make an estimation on the value of a by 
assuming that the potential ([H]) is a solution of differ- 



ential equation (|19p when fi^ 1. Thus, Aq ~ |a and 
consequently a ~ |. By this amount for a, m is negative 
and also it is smaller than n {m— -3.45), as we expected. 
As a result, the following model 



f{R) - i? - nR'' + vIT 



(47) 



can lead to the thawing behavior in the matter dominated 
epoch and late times. Note that this model satisfies the 
condition (|20p . Also, as we required, it's corresponding 
potential in the Einstein frame has a maximum. 



DISCUSSION 

We have considered BD scalar tensor theory in the Ein- 
stein frame. In this frame, BD theory behaves like an in- 
teracting quintessence. It is necessary that the potential 
V{i~p) has a maximum in the region where the scalar field 
rolls in order to have thawing behavior. Also the poten- 
tial should satisfy the condition F < 1. The thawing con- 
dition shows that for non-interacting quintessence 
model, potential should satisfy the condition F w 1 ^19] . 

In the last section, by setting the BD coupling constant 
u) to zero, we have studied the thawing behavior of f{R) 
gravity models in the Einstein frame. It is important 
to note that for power law f{R) gravity models, such 
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as (|45|) . the equation of state parameter of dark energy 
(in the Einstein frame) firstly increases with time and 
then decreases. This behavior is due to the form of the 
corresponding potential of these models in the Einstein 
frame. For these models, A is constant and as we have 
mentioned in section II, co^ evolve as Fig. [TJ The sign 
of the power of R depends on the present day value of 
Lj^p and it's magnitude depends on Aq. As it is clear from 
Fig. [l] choosing different values for Aq leads to different 
values of uj^ at the present day. 

As an example for thawing f{R) models, we have pro- 
posed the model given by (|47[) . The corresponding poten- 
tial has a maximum and the condition F < 1 is satisfied. 
So, in the beginning of the matter dominated era A is 
not approximately constant, (see equation (|34p ). and uj^p 
increases slowly as it has been shown in Fig. |H Also 
this model leads to a nearly flat potential in the late 
times which is satisfactory. At sufficiently late times, the 
interaction term in equation p7|) is negligible and we ex- 
pect our model behaves like a non-interacting thawing 
quintessence [lit . 
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